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We study the Cauchy horizon (CH) singularity of a spherical charged black hole perturbed nonlin- 
early by a self-gravitating massless scalar field. We show numerically that the singularity is weak 
both at the early and at the late sections of the CH, where the focusing of the area coordinate r is 
strong. In the early section the metric is almost Reissner-Nordstrom, and the fields behave accord- 
ing to perturbation analysis. We find exact analytical expressions for the gradients of r and of the 
scalar field, which are valid at both sections. We then verify these analytical results numerically. 
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PACS number(s): 04.70.Bw, 04.20.Dw, 04.25.Dm 

One of the long-standing interesting predictions of 
General Relativity (GR) is the occurrence of spacetime 
singularities inside black holes (BHs). This issue is in- 
triguing, because the laws of physics we currently under- 
stand (e.g. classical GR) are not valid at singularities, 
but some other, as yet unknown, laws take over from 
classical GR and control the structure of singularities. 
Thus, despite many efforts in the last few decades, the 
nature of the spacetime singularities in a generic gravi- 
tational collapse — and, more generally, the final outcome 
of the collapse — are still open questions. 

Until recently, the only known generic singularity was 
the BKL singularity — an oscillatory spacelike singu- 
larity. In the last few years, however, evidence has been 
steadily accumulating that another type of singularity 
forms at the Cauchy horizon (CH) of spinning or charged 
BHs. The features of this new singularity differ dras- 
tically from those of the previously-known singularities 
like, e.g., Schwarzschild or BKL: First, the CH singu- 
larity is null rather than spacelike dj-l); Second, it is 
weak H,^. Namely, the tidal distortion experienced by 
an infalling extended test body is finite (and, moreover, 
is typically negligibly small) as it hits the singularity [|j . 
Yet, curvature scalars diverge there j^,^ (in the spherical 
charged case, this is expressed by mass inflation Q). 

For uncharged spinning BHs (the more realistic case), 
the evidence in favor of this new picture emerges pri- 
marily from a systematic linear and non-linear perturba- 
tive analysis |^J|]. In addition, the local existence and 
genericity of a null weak singularity in solutions of the 
vacuum Einstein equations was demonstrated in Ref. . 
(The compliance of null weak singularities with the con- 
straint equations was demonstrated in Q.) In the case 
of a spherical charged BH, the weakness of the singu- 
larity was first demonstrated in [^. More recently, an 
approximate leading-order analysis |^ confirmed the lo- 
cal consistency of this new picture. 

Despite these recent advances, our understanding of 
the null weak CH singularity is still far from being com- 
plete. In particular, it is important to verify this new pic- 
ture by performing independent, non-perturbative, anal- 



yses. This motivates one to employ numerical tools to 
study the structure of the CH singularity. The numer- 
ical simulation of spinning BHs is difficult, as they are 
non-spherical. One is thus led to study, numerically, the 
inner structure of a spherical charged BH; hopefully, it 
may serve as a useful toy model for a spinning BH. 

The first numerical analysis of a perturbed charged 
BH's interior was carried out by Gnedin and Gnedin 
|]l0| , who analyzed the spherically-symmetric gravita- 
tional collapse of a self-gravitating scalar field over a 
charged background. The coordinates (and numerical 
grid) used there, however, do not allow getting even close 
to the CH. More recently, Brady and Smith (BS) |l|] 
numerically explored the mass-inflation singularity in- 
side a spherical charged BH perturbed nonlinearly by a 
scalar field. This analysis confirmed several aspects of 
the above new picture: It demonstrated the existence of 
a null singularity at the CH, where the mass function 
TO diverges but the radial Schwarzschild coordenate r is 
nonzero. The quantity r was found to decrease mono- 
tonically with increasing retarded time along the CH, 
due to the nonlinear focusing, until it shrinks to zero (at 
which point the singularity becomes spacelike). It also 
provided evidence for the weakness of the singularity. De- 
spite its remarkable achievements, however, this analysis 
left one important issue unresolved: To what extent is the 
perturbative approach applicable at (and near) the CH 
singularity? BS reported on an inconsistency with the 
predictions of perturbation analysis, manifested by the 
non-zero value of a (see pT[|), namely a finite deviation 
of the power-law indices from the integer values predicted 
by perturbative analyses. This issue is crucial, because 
for reahstic (i.e., spinning and uncharged) black holes the 
only direct evidence at present for the actual occurrence 
of a null weak singularity stems from the perturbative 
analysis A failure of the perturbative approach in 
the spherical charged case would therefore cast doubts 
on our understanding of realistic black holes' interiors. 

In this paper we report on a numerical and analytical 
investigation, which we carried out in order to answer 
this and other questions. We consider the model of a 
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spherical charged black hole non-linearly perturbed by 
a spherically-symmetric, self-gravitating, neutral, mass- 
less scalar field $ (the same model as in | |To| , [ll[ | ) . We 
shall first present our numerical results which show that 
the CH singularity is indeed weak — not only in the early 
part of the CH, but also all the way down to the point of 
complete focusing (r = 0). Then, we study the asymp- 
totic behavior of perturbations at the early part of the 
CH, and demonstrate the full compliance with the predic- 
tions of perturbation theory (in particular, we find that 
cr = 0). In addition, we shall show that despite the diver- 
gence of m, the metric functions at the asymptotically- 
early part of the CH are remarkably close to the un- 
perturbed Reissner-Nordstrom (RN) metric functions — 
which is again a prediction of the perturbation analysis 
(according to the latter, the metric perturbations should 
vanish at the asymptotic past "edge" of the CH, despite 
the divergence of the curvature [D). Then, we shall an- 
alyze the behavior of the blue-shift factors and <I>_^, 
as a function of u, at the late (i.e. strong-focusing) part 
of the CH, where u and v are ingoing and outgoing null 
coordinates, respectively (see below). The asymptotic 
behavior of r^y and <i>^^, is essential, because it is primar- 
ily these entities that are responsible to the divergence 
of curvature at the CH. We shall present exact analytic 
expressions for these entities, and verify them numeri- 
cally. The expression we obtain for in addition to 
our numerical results, shows that a vanishes not only at 
the early part of the CH, but also everywhere along it. 

We write the general spherically-symmetric line ele- 
ment in double-null coordinates. 



Spacelike 
singularity 



-f{u, v) du dv + r'^{u, v) d^l^ , 



(1) 



where dil^ is the unit two-sphere. As the source term 
for the Einstein equations we take the contributions of 
both the scalar field and the (sourceless) spherically- 
symmetric electric field (see for more details). The 
dynamical field equations are 
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where the constant Q is the electric charge. Equations 
(h)~(q) ^-re supplemented by two constraint equations: 



r,uu - (ln/),„r„ -|- r{^^uf = 



r,vv - (ln/)^t,r 1, -|- r($,^)^ = 0. 



(5) 



(6) 




Initial pulse 



FIG. 1. The Penrose diagram of the simulated spacetime. 
Singularities are displayed by thick lines, and the character- 
istic hypersurface on which the initial data are specified are 
marked with dashed lines. The various fields are probed along 
outgoing (such as 1 — in the early sections — and 2 — in the late 
sections) and ingoing (such as 3) null rays. It is still unknown 
whether there is a continuation of the spacetime manifold be- 
yond the CH. Our domain of integration cannot include the 
entire spacelike portion of the singularity. These uncovered 
areas are marked by dashed thick lines. 

It turns out that it is advantageous to substitute 
f{u, v) = 2e^*'-"'^^ for the numerical integration near the 
CH. Our initial- value setup is described in Ref . : The 
geometry is initially RN, with initial mass Mq = 1 and 
charge Q, and no scalar field. At some moment u, how- 
ever, it is modified by an ingoing scalar-field pulse of a 
squared-sine shape with amplitude A. $ vanishes ev- 
erywhere on the initial surface except in a finite range 
vi < V < V2- The results presented below relate to 
vi = 10, V2 = 20, Q = 0.95 and A = 8 X 10"^ un- 
less stated otherwise. In this case, due to the scalar- 
field energy, the BH's external mass approaches the fi- 
nal mass Mf « 1.4. (We also checked other values of 
0.5 < Q/Mf < .99 and A, and obtained similar results.) 
Note that our outgoing initial null hypersurface is lo- 
cated outside the event horizon (EH) (unlike in Ref. |0 ) . 
Therefore, we do not have to make any assumption about 
the inverse power tails at the EH; these are created au- 
tomatically by the dynamical evolution. Our numerical 
scheme is essentially the same as described in Ref. |lj] 
(there are few modifications, which will be described else- 
where [|lj]): It is based on free evolution in double- null 
coordinates. The code is stable and second-order accu- 
rate p2[ . Our numerical setup is displayed in Figure 1, 
imbedded in the Penrose diagram of the simulated space- 
time. 

The null coordinates u and v are defined in Ref. p2[ : 
They are taken to be linear with r on the two character- 
istic initial segments. For the presentation and interpre- 
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tation of the results, we shall occasionally use other types 
of double- null coordinates: Eddington-like coordinates Ue 
and Ve, and Kruskal-like coordinates (of the inner hori- 
zon of RN) U and V. (In the perturbed spacetime, the 
Eddington-like and Kruskal-like coordinates are defined 
with respect to the asymptotic "vertex" region, where 



the metric perturbations vanish asymptotically; see |14|) 
Recall that at w » M/, v is closely related to Ve 

Our numerical simulations confirm the presence of a 
null singularity at the CH, where m diverges and r is 
nonzero. Along the CH singularity, r decreases mono- 
tonically, until it shrinks to zero, at which point the sin- 
gularity becomes spacelike. This situation was already 
found numerically by BS 
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FIG. 2. Metric functions in the early sections of the CH. 
g [(A) and (B)] and r [(C) and (D)] as functions of Iog(f/V) 
along an outgoing ray and an ingoing ray ai v — 80, respec- 
tively. The solid lines are the RN values (computed according 
to the BH's external parameters at « = 75), and the circles 
represent the numerical values. (E): The mass function along 
the same outgoing ray as (A) and (C). Here we took A^ = 80 
and A = 5 X 10^^, which corresponds to Mf ~ 1.1. 

We shall first discuss the early part of the CH singular- 
ity, i.e., the part where the focusing of r is still negligible. 
Our first goal is to demonstrate that the singularity is 
weak. In terms of the double- null metric (^), the singu- 
larity will be weak if coordinates u{u), v{v) can be chosen 
such that both r and are finite and nonzero at the CH. 
The numerical analysis by BS already demonstrated the 
finiteness of r, which we recover in our results. (Note that 
r is independent of the choice of the null coordinates.) 
Figure 2(A) displays the metric function g = —2guv in 
Kruskal-like coordinates U, V along an outgoing null ray 
that intersects the early section of the CH singularity. 
The CH is located at V — (corresponding to w — > oo). 
This figure clearly demonstrates the finiteness of guV: 
from which the weakness of the singularity follows. 

The perturbation analysis also predicts that both the 
scalar field and the metric perturbations will be arbitrar- 



ily small at the early section of the CH. In other words, 
both metric functions r and g should be arbitrarily close 
to the corresponding RN metric functions. This behavior 
is indeed demonstrated in Fig. 1. In this figure, g and r 
are displayed along lines u = const [Figs. 2(A) and 2(B)] 
and V = const [Figs. 2(C) and 2(D)]. The similarity of 
the analytic RN functions and the numerically-obtained 
functions of the perturbed spacetime is remarkable. (We 
emphasize, though, that despite the similarity in the val- 
ues of the metric functions to RN, our geometry is dras- 
tically different from that of RN, in the sense that in our 
case curvature blows up at the CH.) 
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FIG. 3. Fields in the early sections of the CH. (A): * 
(solid), ^l/.i, (dashed), and r^^ (dash-dotted) along an outgoing 
ray, as a function of v. (B): R, and (C): log m, along the same 
outgoing ray, as a function of v. (D): \l/ along an ingoing ray, 
as a function of U. 

The validity of the perturbation analysis is checked in 
more detail in Fig. 3, which displays the power-law be- 
havior of 5* = r<I>, 5*^1,, and r^^, near the early section of 
the CH. Fig. 3(A) displays these entities along an out- 
going ray. Both the ringing and the power-law tails are 
apparent (the graph of ^' approaches a nonzero limiting 
value, due to the finite value of U on that ray). Fig. 3(B) 
displays R = where = «'™[2(M//g)2 - 1] 

is the asymptotic form of y as predicted by linear per- 
turbation theory |l^, and is the (w-dependent) value 
of at the EH. Numerically, we find that at large v, 
R tends asymptotically to 1 — it deviates from unity by 
no more than 2 x 10"'^ at v = 600. A similar result was 
obtained for various values of Q and A. We thus find 
no evidence for charge-dependent internal power-law in- 
dices, and we conclude that the parameter a of Ref. |Tl| ] 
vanishes identically. The behavior of m along the same 
outgoing ray is shown in Fig. 3(C). Clearly, the exponen- 
tial growth of m does not affect the validity of the per- 
turbation analysis. Finally, Fig. 3(D) displays 'if along 
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an ingoing null ray at very large v. The local power index 
(see [H]) of Uf, is found numerically to be 3.1, in a re- 
markable agreement with the value 3 predicted by linear 
perturbation theory. We conclude that the linear pertur- 
bation analysis describes the evolution of perturbation 
near the vertex well, despite the divergence of m (and 
the curvature) . This confirms the prediction of nonlinear 
perturbation analysis Q|, namely that nonlinear pertur- 
bations will be negligible at the early part of the CH, 
compared to the linear perturbations. 

We turn now to explore the late part of the CH sin- 
gularity, where focusing is strong. First, we numerically 
verify the weakness of the singularity in this part too. 
Figure 4(A) shows G = —2guv along an outgoing null 
ray in the late part of the CH, where the value of r has 
shrunk to 10% of the value it had when the CH was first 
formed. We present here the results for various values of 
the grid-parameter N (see @), in order to demonstrate 
the second-order numerical convergence, guv approaches 
a finite value at the CH {v — > oo). At the same time, the 
mass function (and curvature) grows exponentially with 
V [see Fig. 4(B)]. We conclude that the entire null CH 
singularity is weak, even at the region of strong focusing. 



[2(My/Q)-l] vE- 



EH 



(8) 





Here, k- is the surface gravity at the RN inner horizon 
with parameters M/ and Q. For convenience, Eqs. (0) 
and (||) are expressed in terms of 5*. Note that (^) is 
invariant to a gauge transformation v — > v{v), whereas 
(0) is not. [Eq. (0) refers explicitly to the derivatives 
with respect to the Eddington-like coordinate Vg.] 

In order to verify this prediction, we calculated nu- 
merically R = ^',„/{^'™[2(M//Q)2 - 1]} and P = 
— (r^) ^ / [(1/k_)^'_„^^] , as functions of along an out- 
going ray located at a region of 90% focusing of the CH. 
The results, presented in Fig. 4(C), are in excellent agree- 
ment with the above theoretical prediction, R ^ 1 ~ P. 

To summarize, we have confirmed, numerically (for 
a spherical charged black hole) the main predictions of 
the perturbation analysis (which apply both to spherical 
charged and nonspherical, spinning BHs): The null sin- 
gularity at the CH is found to be weak. In the asymptotic 
early section of the CH the metric functions approach 
arbitrarily close to the corresponding metric functions in 
RN, and, moreover, the perturbations are well described 
by the linear perturbation analysis. This confirms the 
conclusions of the perturbation analysis 0| , that the non- 
linear effects are negligible at the early section of the CH 
singularity. In addition, we analytically derived exact ex- 
pressions for the diverging blue-shift factors r.„ and ^^v, 
which are valid everywhere along the contracting CH. 

I am indebted to Amos Ori for numerous invaluable 
discussions and useful comments. 



FIG. 4. The nonlinear regime (contraction of the CH of 
90%): (A) G along an outgoing null ray, as a function of v. 
Dotted: iV = 20, dash-dotted: = 40, dashed: N = 80, and 
solid line: N = 160. (B) log m along the same outgoing ray. 
(C): R (dotted) and P (solid) along the same outgoing null 
ray, as a function of v. 

Next we study, analytically, the behavior of the blue- 
shift factors and <I>^^ along the contracting CH. Here, 
we shall present the results; the full derivation will be 
presented in Ref. [Q. The field equations (||) and (E|) 
can be integrated exactly along the CH singularity [[l4[ . 
For r_y we find 



and for 'i! ^ we find 
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